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The transport of sputtered aluminum inside a multi frequency capacitively coupled plasma chamber 
is simulated by means of a kinetic test multi-particle approach. A novel consistent set of scattering 
parameters obtained for a modified variable hard sphere collision model is presented for both argon 
and aluminum. An angular dependent Thompson energy distribution is fitted to results from Monte 
Carlo simulations and used for the kinetic simulation of the transport of sputtered aluminum. For 
the proposed configuration the transport of sputtered particles is characterized under typical process 
conditions at a gas pressure of p = 0.5 Pa. It is found that - due to the peculiar geometric conditions 
- the transport can be understood in a one dimensional picture, governed by the interaction of the 
imposed and backscattered particle fluxes. It is shown that the precise geometric features play an 
important role only in proximity to the electrode edges, where the effect of backscattering from the 
outside chamber volume becomes the governing mechanism. 


I. INTRODUCTION 

Capacitively coupled plasmas are commonly used for 
plasma enhanced chemical vapor deposition (PECVD) 
applications [T] . As recently demonstrated by Bienholz et 
al. [2]-[5] capacitive discharges operated at multiple driv¬ 
ing frequencies can be also employed for physical vapor 
deposition (PVD) and sputtering applications for the de¬ 
position of metallic, ceramic, and/or magnetic thin films. 
In contrast to cathodic arc evaporation [6], magnetically 
enhanced PVD processes such as the widely used direct 
current magnetron sputtering (dcMS) [7], or the recently 
technologically evolving high power impulse magnetron 
sputtering (HiPIMS or HPPMS) [HIH], multi-frequency 
capacitively coupled plasmas (MFCCP) for sputtering 
applications have the advantage of a diverse choice of ma¬ 
terials, including non-conducting and ferromagnetic tar¬ 
get materials [2]. For these materials DC driving sources 
and/or magnetic field enhancement may not be a suit¬ 
able option. Despite the fact that the MFCCP reactor 
proposed by Bienholz et al. has been experimentally char¬ 
acterized in great detail, so far only a preliminary expla¬ 
nation of the transport phenomena of sputtered particles 
is provided [3] . For the interpretation of some of the ex¬ 
perimental findings there is a lack of reliable simulation 
data. 

To close the gap between experiments and theory the 
transport of sputtered particles is studied under realis¬ 
tic discharge conditions. In order to characterize the 
geometrical constraints for the deposition behavior the 
transport of sputtered aluminum is first simulated within 
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a fully three dimensional geometry. The interaction of 
the sputtered aluminum with both the background ar¬ 
gon gas and the backscattered aluminum is subsequently 
analyzed using a reduced one dimensional setup. The 
comparison of the two models reveals the suitability and 
limitations of the reduced ID model governing the trans¬ 
port and thus the density and mean velocity of the sput¬ 
tered particles. 


For the proposed theoretical approach a strongly modi¬ 
fied version of the direct simulation Monte Carlo (DSMC) 
code dsmcFoam [9] is used. (dsmcFoam is part of the 
open source software package OpenFOAM [lO]). Al¬ 
though the code in principle allows to self-consistently 
simulate the background flow field as well as the sput¬ 
tered particle transport, in this work, for the sake of 
computational feasibility, only the transport of sputtered 
particles is simulated using a variant of the test particle 
method (TPM). This is necessary because the sputtered 
aluminum particles constitute a trace minority in the 
overall gas composition, mainly the argon background. 
The TPM has been used for numerous studies of sputter¬ 
ing processes as reported by a number of researchers and 
research groups [IIH23]. Several different implementa¬ 
tions have been referenced in the literature, for example 
a single particle approach m, or a treatment of inde¬ 
pendent particles in a parallel fashion following the test 
multi-particle method (TMPM) [151 [16]. either case, 
three aspects are of particular importance: i) the details 
of the numerical implementation of the TPM itself, ii) 
the underlying collision model for the particle interac¬ 
tion, and iii) the angular and energy distribution of ini¬ 
tially injected sputtered particles. All of these aspe cts a re 
innovatively reconsidered in this work (cf., section III). 

The manuscript is organized as follows: After a brief 
description of the (experimental) discharge configuration 
and the operating conditions in section [II] a review of the 


employed numerical scheme is presented in section III 


Particular emphasis is placed on the collisional interac¬ 
tion of sputtered particles with the background gas as 
well as on their initial energy distribution after sputter- 
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Figure 1. Sectional drawing of the MFC CP vacuum chamber. 
The driven top electrode is the target, while a substrate can 
be mounted at the opposite bottom electrode. The latter as 
well as the chamber walls are grounded. The red dashed frame 
indicates the electrode gap region. 

ing occurs. Then in section |IV[ firstly, simulation re¬ 
sults for the realistic three dimensional MFC CP vacuum 
chamber are discussed with the focus on the geometric 
constraints affecting the details of the particle transport. 
Secondly, simulation results from a reduced one dimen¬ 
sional model are discussed highlighting the fundamental 
transport mechanism. A parameter study of the den¬ 
sity and the mean velocity profiles for the sputtered alu¬ 
minum is presented. The nature of the governing trans¬ 
port process is discussed on the basis of the particles’ ve¬ 
locity distribution functions (VDF). Finally, in section [V| 
the results are summarized and conclusions are drawn. 


II. DISCHARGE SETUP 

To obtain a principle understanding of the subject under 
investigation in this work, a brief review of the discharge 
setup previously described and investigated by Bienholz 
et al. is instructive 1315]. A sectional drawing of the 
vacuum chamber is depicted in figure The stainless 
steel vessel forms a cuboid of 80 cm side lengths and a 
height of 65 cm. Two opposing circular electrodes of di¬ 
ameter VF = 50 cm are separated by a gap of distance 
D = 7.5 cm. While the upper (target) electrode is driven 
with multiple radio-frequencies (RF), the lower (sub¬ 
strate) electrode, as well as the remaining parts of the 
reactor chamber (i.e., the walls) are electrically grounded. 
Consequently, under operation this geometrically asym¬ 
metric discharge exhibits a substantial self-bias voltage. 
The discharge is typically simultaneously driven by three 
independent power supplies at frequencies of 13.56, 27.12 
and 60 MHz. Taking advantage of the electrical asymme¬ 
try effect (EAE) [24l [25] the MECCP has been demon¬ 
strated to be an ideal candidate as a capacitively cou¬ 
pled sputtering system. The target electrode is made of 


aluminum and is subject to undergo sputtering by ion 
bombardment. The remaining components are made of 
stainless steel. Since there is no necessity to magnetically 
enhance the plasma, besides the inherent plasma profile 
over the electrode no additional lateral inhomogeneity 
is introduced due to magnetic fields (e.g., a race-track). 
Consequently, advantage can be taken from a nearly ho¬ 
mogeneous target utilization over a quite large area of 
Atarget ~ 1963 CITL^. 

Bienholz et al. provide a sophisticated experimental 
characterization based on a variety of diagnostic tech¬ 
niques (i.e., optical emission spectroscopy, Langmuir 
probe, and micro balance deposition measurements). 
They further specify the parameters for a reference sce¬ 
nario to be investigated in this work: In steady state 
operation the gas temperature T = 650 K is measured 
with a confidence interval of ±50 K. At a gas pres¬ 
sure of p = 0.5 Pa the background gas density (ar¬ 
gon) then is nAr = 5.572 x 10^^ cm“^. The energy of 
monoenergetic argon ions impinging the target, Ei, is 
taken from the measured self-bias voltage Vsb = 437 V. 
The assumption of monoenergetic ions is justified be¬ 
cause the ion plasma frequency cjpi « 28 x 10^ s“^ is 
significantly smaller than the lowest driving frequency 
<^rf, min = 27r X 13.56 MHz ^ 85.2 x 10^ s“^. Ions there¬ 
fore experience only the average electric field within the 
boundary sheath in front of the electrodes. Consequently 
they are not modulated. Moreover, under the assump¬ 
tion of a low collisionality (which is the case here), the 
ion energy distribution is not substantially broadened. 

Advantage can be taken from simulations of the sput¬ 
tering process itself using the code TRIDYN [261 EZj. 
The angular and energy distributions of sputtered par¬ 
ticles regarded in this work have been obtained by Bi¬ 
enholz et al. like so |3|. The flux of aluminum particles 
sputtered from the target is estimated from the exper¬ 
iment as Fo = 1.6 X 10^^cm“^s“^ [5]. A homogeneous 
flux of sputtered particles over the target surface can 
be assumed, because of the quite homogeneous plasma 
density in front. Moreover, based on retarding field mea¬ 
surements performed by Ries et al. [28] resputtering from 
the substrate surface can be neglected. From the low ion 
energies of Einc ^ 50 eV at the substrate the sputtering 
yield is estimated to be less than a few percent. Finally, 
at all surfaces complete adsorption of the impinging alu¬ 
minum can be assumed (with a sticking coefficient 5 = 1). 
Although rather simple, this is a good approximation for 
(low) energetic aluminum atoms (< 25 eV), regardless of 
their angular distribution [29] . 

Within the volume typical aluminum densities of ap¬ 
proximately 10^^ cm“^ result. In comparison, the back¬ 
ground density is on the order of 5 x 10^^ cm“^. From 
the ratio of the two (i.e., 2 x 10“^) it can be reasoned that 
aluminum is merely contained in traces. In consequence, 
even though the average energy of sputtered particles is 
roughly two orders of magnitude higher than the back¬ 
ground’s thermal energy (k^T « 56 meV), gas heating 
and thus rarefaction of the background gas remain negli- 
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gible. Effects like a sputtering wind m can therefore be 
excluded. (Note that the contrary is typically the case 
for HiPIMS. In that case a substantial rarefaction of the 
background may occur due to high metal densities com¬ 
parable to the density of the process gas, paired with a 
high degree of ionization. The proposed model in the 
present form is inappropriate for this situation.) 

In addition to the experimental analysis, Bienholz et 
al. also provide simulation results obtained from a one di¬ 
mensional TPM code based on collisional data of Kuwata 
et al. m- Despite the strong effort, these theoretical 
results do not seem to explain what is observed experi¬ 
mentally. In particular there are some quite contradic¬ 
tory discrepancies in the sputtered particle densities and 
fluxes. The goal of this work is to provide a detailed the¬ 
oretical analysis of the underlying transport mechanism 
in order to fully understand and explain these seemingly 
contradictory results. 


III. NUMERICAL METHOD 


As previously introduced, in this work a strongly modi¬ 
fied version of the open-source DSMC solver dsmcFoam is 
employed 0110 ]. The original version of the package pro¬ 
vides the basic functionality in terms of a Lagrangian par¬ 
ticle motion (with one super-particle representing a given 
number of physical particles), binary collisions based on 
the variable hard sphere (VHS) model, as well as free- 
stream boundary conditions. Its capability for the sim¬ 
ulation of the particle transport in sputtering systems is 
however limited. A number of necessary modifications 
have been applied to the code which comprise of the im¬ 
plementation of the Ml binary collision model [32], inlet 
and wall boundary conditions appropriate to accurately 
mimic sputtering from a given target surface assuming a 
modified Thompson distribution [33], and species selec¬ 
tive particle sticking to the walls. The specifics of these 
modifications are detailed in subsections IIII Al and IIII Bl 
For the transport simulation of trace species a sim¬ 
plified computational procedure can be adopted. The 
minority particles solely interact with a constant, non 
stationary background, but not among themselves. This 
(or an equivalent) assumption is common to test-particle 
methods in general. The background is then completely 
specified by its number density ribg and temperature T. 
The proposed TMPM model differs significantly from 
the previously described TMPM models [la US]. In 
our model a simplified no-time counter (NTC) method 
paired with a constant time-stepping scheme is used. For 
all super-particles (also referred to as simulators) equal 
weights are assumed. Correspondingly, the number of 
physical particles represented by one super-particle is 
w = nVc/Nc, with the trace species number density n, 
the cell volume Vc and the number of super-particles per 
cell Nc. The probability for a trace species simulator i to 
undergo collision with the background can be specified 


as [16] 


^i,bg — '^bg (^Tld)-At, 


( 1 ) 


where id = \v[ — Vhg\ is the relative velocity between 
the simulator i and the background sampled from a 
Maxwellian distribution with temperature T. ax = 
crT(id) is the velocity dependent total collision cross- 
section integrated up to a finite cut-off 6max of the impact 
parameter h (cf., section III A). At is the discrete simu¬ 
lation time step. In the fashion of the NTC method, the 
collision probability can be rewritten as 


-Pi,bg — Acand 


■ 1 ■ 



_max((7TV;)ibg_ 


( 2 ) 


with the number of collision candidates per time step 
given by 


Ac and 


ribgAc max {crTVr)-^^^ At 


(3) 


The first pair of parentheses in eq. (§ represents the 
probability of randomly drawing one out of Nc simula¬ 
tors, the second pair of parentheses is a normalized col¬ 
lision probability for this choice. The computational ef¬ 
ficiency can thus be largely improved by selecting only 
Acand candidates out of the Nc simulators per cell. This 
representation is physically equivalent to the TMPM as 
described in [Tgnni, only that it is significantly more ef¬ 
ficient. The restriction of a maximum allowed time step 
At to ensure the probability of a collision to be Pcoi < 1, 
is not truly circumvented as also Acand < Ac has to be 
ensured. Thus At ^ Tc is a priori chosen to be much 
smaller than the mean collision time Tc [35l [36] . 

On collision, the actual scattering is evaluated 
The pre-collision velocities are first transformed from the 
laboratory into the co-moving center of mass frame. In 
the co-moving frame local reference coordinates are then 
defined. Next, the impact parameter b = 6niax\/^ 
is selected (uniformly distributed over the cross-section) 
with a uniformly distributed random number P/ G [0,1] 
[nnnEz]. Along with the impact parameter also the 
scattering angle x(6) is specified using in our case the 
Ml model (to be discussed next). The azimuthal angle 
-0 = 27rRf is drawn (with another independent random 
number Rf) to completely specify the post-collision scat¬ 
tering direction in the local reference frame. The post¬ 
collision velocities are finally obtained in the laboratory 
frame (taking into account the rotation and change in 
coordinates). 


A. Ml Collision Parameters 

In favor of the ordinary VHS model which comes with 
dsmcFoam the Ml collision model proposed by Morokoff 
and Kersch is used [32]. (Therein an elaborated discus¬ 
sion of the Ml model is provided.) As reasoned, com¬ 
pared to the VHS model it gives a more accurate de¬ 
scription of the collision dynamics at the cost of a slightly 
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Table L Ml collision parameters obtained for argon and alu¬ 
minum. 




Argon 

Aluminum 

dref,i 

(A) 

4.614 

4.151 

LUi 

(-) 

0.721 

0.72 


more sophisticated evaluation. This comes from the scat¬ 
tering direction being not isotropic on a sphere anymore. 
Despite that the Ml collision model can be used with es¬ 
sentially the same parameters as the VHS model (scaled 
by a factor of a/ 4/3 [15] ) , in this work a novel set of colli- 
sional parameters for all molecular interactions included 
in the simulation is proposed. This is necessary because 
for the collision processes of interest so far only rather 
inconsistent VHS or Ml data sets are available (in par¬ 
ticular for metal species). In order to obtain a consistent 
set of cross-sectional data, the scattering (or deflection) 
angle |37| 



(4) 


is first calculated from an analytic interaction poten¬ 
tial 4>(r) as a function of the impact parameter b and 
the relative velocity V- (writing the collisional energy as 
^{Vr) = ^nirV^). The distance of closest approach is va 
and the reduced mass is calculated from the individual 
particle masses = 17111712 /{mi 4 - 7712 ). The interac¬ 
tion potential is 4>(r) = Aex.p {—Br) based on a Born- 
Mayer form [38] . The well-documented set of parameters 
of Abrahamson is used [39] . The scattering dependence is 
obtained for energies 5 (V) = {0.01, 0.1, 1, 10, 100} eV. 
Over the range of these energies the total deviation is 
minimized using the method of least squares. Corre¬ 
spondingly, the Ml model’s linear relationship 

x{b,Vr)=7r[l-b/d,{Vr)] (5) 

is used. The cut-off impact parameter 6rnax equals di(Vr), 
the velocity dependent diameter of species i. Following 
the commonly used temperature power law for the vis¬ 
cosity /i(T) oc (T/Tref)^' [34l|35j|^ the latter reads 


di{Vr^ — <^ref,i 


/ 






1 

r(5/2-u;i)' 


(6) 


Figure shows the obtained scattering angle 
x{b, e{Vr)\ diS well as the linear approximation for ar¬ 
gon (Fig.j^) and aluminum (Fig. |^). The obtained Ml 
parameters drei,i and cji are listed in table Clearly, the 
linear fits resemble the interaction kinetics of the much 
more complex interaction potential reasonably well. For 
all energies the overall functional shape is similar. A 
systematic deviation is observed for impact parameters 



b(A) 



b(A) 

Figure 2. Scattering angle y plotted as a function of the im¬ 
pact parameter b for the linear dependence of the Ml model 
(solid lines) and the analytic calculations (dotted lines) for 
different interaction energies, a) Argon (above) and b) alu¬ 
minum (below). 


b > di{Vr) due to the assumption of a finite distance of 
molecular interaction. For b < di{Vr) the dependence of 
the scattering angle on the interaction energy is captured 
appropriately over a wide range of energies. Only for very 
low (5 < 10 meV) and quite high {e > 100 eV) energies 
the Ml model more substantially deviates from the ana¬ 
lytic calculations. This is a limitation of the power law 
assumption. Depending on the energy range of interest 
a more appropriate energy dependence may be chosen, 
possibly involving a more focused consideration of the 
kinetics of molecular interaction. 

Of major importance for this work are combined col¬ 
lisions between argon and aluminum. For the ana¬ 
lytic interaction potential as well as for the Ml model 
the corresponding collision parameters can be evaluated 
based on the parameters specified for the individual 
species. For the analyt ic interaction potential a com¬ 
bining rule 4>(r) = Y^4>Ar(^)^Ai(^) is used [39l [40] . 
For the Ml model a linear interpolation dref,Ar-Ai = 
(<^ref,Ar + (iref,Ai)/2 and c^Ar-Ai = (^Ar + cjai)/2 is per¬ 
formed [ 34 ]. The corresponding results obtained from 
the analytic interaction potential as well as the Ml de¬ 
pendence for Ar-Al collisions are displayed in Fig. |3}i. It 
is important to notice that also in the combined case the 
scattering characteristics are captured well - within the 
limitations of the individual argon/aluminum fits. 

The differential cross-section for an Ar-Al collision as 
calculated from the scattering dependence of the Ml 
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Figure 3. a) Scattering angle x of Ar-Al collisions plotted as 
a function of the impact parameter b for the linear depen¬ 
dence of the Ml model (solid lines) and the corresponding 
analytic calculations (dotted lines) for different interaction 
energies, b) The corresponding Ar-Al differential collision 
cross-sections based on eq. 0- 


model im 


o'(x) 


b db 


sinx 


dx 


dHVr) q _ x\ 

TT sin X ^ TT / 


(7) 


is shown in Fig. [^. Its angular dependence forthrightly 
reflects the anisotropy of the underlying collision model 
(specifically, there is no symmetry about 7r/2). For glanc¬ 
ing collisions (i.e., b <ii(K) and thus X ^ 0) the dif¬ 
ferential collision cross-section diverges. This stems from 
sin X in the denominator and reflects the vanishing of the 
differential dft = sinxdxdij. The singularity at x = 
is canceled. Notably, integrals over a{x)dQ remain fi¬ 
nite. In consequence the total collision cross-section ctt, 
the momentum (or viscosity) cross-section ctm, and the 
diffusion cross-section ctd can be obtained [351EZ]. The 
total collision cross-section which enters into eqs. 0 to 
^ finally reads ctt = 7 rdf{Vr). 


B. Sputtered Particle Distribution Function 

The angular and energy distribution of particles sput¬ 
tered of various kinds of surfaces has undergone exten¬ 
sive analyses in the past decades starting from the work 
of Thompson and Sigmund jiaiia] and including numer¬ 
ous other works mEsmiHii]. In this paper a rather 
empirical fitting approach is pursued where the angular 
and energy dependences of the respective distributions 
are compared and fitted to results reported by Bienholz 
et al. obtained from TRIDYN simulations [3]. Instead of 
directly using their Monte Carlo results, analytic expres¬ 
sions are used in the TMPM simulations to initially spec¬ 
ify the sputtered particle velocities. For this purpose, the 



Figure 4. Angular distribution of sputtered aluminum ob¬ 
tained from TRIDYN simulations (solid black line) by Bien¬ 
holz et al. [3] and analytic fit using a cosine distribution eq. 0 
(dashed red line). 


normalized fiux distribution of sputtered particles may be 
denoted by [22] 

f{E,0,(p)dEdn = f{E\0)dE x f{0) sinOdO x f{ip)dip, 

(8) 

with the solid angle dCl = sin OdOdcp. (p is the uni¬ 
formly distributed azimuthal emission angle with f{(p) = 
(27r)“^. The energy distribution f{E\0) is assumed to be 
dependent on the polar angle of emission 0 which itself 
is measured from the surface normal (0 < 0 < 7r/2). For 
the latter a cosine-law [33l |48] is assumed with 

f{0) =5cos^'((9). (9) 


5 is a normalization factor such that f{ 0 ) sin OdO = 
1. > 0 is a parameter shaping the distribution to 

be over or under-cosine, respectively. In order to cap¬ 
ture the angular distribution obtained by Bienholz et al. 
a = 1.1111 was fitted making the distribution slightly 
over-cosine mi ssi. The corresponding distribution is 
displayed in figure The excellent agreement of eq. ^ 
with the Monte Carlo data is evident. 

As concerns the particles’ energy distribution often the 
angle-independent Thompson distribution is used [42] . It 
is, however, only a rough estimate as observed in exper¬ 
iments and simulations mus]. Based on the work of 
Stepanova and Dew |33| an angle-resolved energy distri¬ 
bution 


f{E\e) = c{e) 


E 

{E + f7)3-2™ 


X exp 




uee 

U + AF^inc 


rmi UEEsliO) 

yTTT/t -E'inc 


n- 

Q' 


( 10 ) 


can be utilized - that is, the conditional probability given 
0. The angular dependence is specified by the shape 
function S 0 {O) = cos^, which is clearly asymmetric with 
respect to 7r/4. The given angular dependence strictly 
only holds for ions under normal incidence. The case of 
ions under oblique incidence is a different matter jH] |47| . 
This topic is beyond the scope of this work, as the as¬ 
sumption of normal incidence is well justified due to a 
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Figure 5. Energy distribution of sputtered aluminum ob¬ 
tained via TRIDYN simulations (solid black lines) by Bien- 
holz et al. [3] and analytic fit using a modified Thompson 
distribution eq. ( |10| ) (dashed red lines) for different emission 
angles 0. 


high sheath potential and a low collisionality for the sit¬ 
uation of interest. The parameters specifying the dis¬ 
tribution of eq. (10) can be divided into two groups: i) 
the rather empirical parameters m = 0.212, A = 0.14, 
n = 0.5, A = 13, g' = 2 — mt/(4mi), and Q = 0.55 
which govern the functional shape of the distribution. 
These parameters are specified in [33] and/or fitted to 
the TRIDYN results, ii) A set of parameters which are 
physically justified. That is, = E^b,Ai = 3.36 eV is 
commonly chosen as the binding energy, the energy of 
ions incident on the target is E^inc = eRsB, the incident 
ion mass is rrii = niAr and the target atomic mass is 
rrit = mAi- C{0) is again used for normalization and 
chosen such that f(El0)dE = 1. 

The expression referenced for the angle-dependent en¬ 
ergy distribution /(Ejd) seems to be appropriate for most 
sputtering applications [22j ESI 113 . However, when com¬ 
paring the TRIDYN results to the proposed expression, 
the angular dependence is captured well only for emis¬ 
sion angles d < 7r/4. Above 0 ~ 7r/4, for the Monte 
Carlo results the high energy tail of the distribution 
drops, while the cos 6>-dependence in the original expres¬ 
sion yields a monotonic rise of the latter. To better re¬ 
flect the Monte Carlo data, instead, a shape function 
S 0 (d) = 1 — /csin(2^) is suggested. Its symmetry about 
7r/4 allows to better capture the observed angular de¬ 
pendence, where k = 0.5 determines its strength. It is 
important to note that due to the choice of the exponents 
q and Q the symmetry about 7r/4 is broken, despite the 
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Figure 6. Sectional drawing of the reactor chamber zoomed to 
the electrode gap region indicated by the dashed red frame in 
figure The aluminum density n is displayed. The dashed 
black lines indicate the positions to be considered later in 
figure 


symmetric choice of se{0). The angle-resolved energy dis¬ 
tribution f{E\0) together with the Monte Carlo results of 
Bienholz et al. is depicted in figure for emission angles 
0 = {0,7r/4,7r/2}. As most evident from the high energy 
tail, for emission angles larger than 0 > it/A the distribu¬ 
tion drops again. This is well captured by the modified 
shape function. Even for the improved angular depen¬ 
dence some deviations are observed. In particular for 
energies from a few up to < 30 eV the analytic expres¬ 
sion slightly overestimates the Monte Carlo results. The 
overall trend however is captured much better, compared 
to a monotonically rising high energy tail assumed for the 
original shape function or the classical angle-independent 
Thompson distribution. Without going into detail it can 
be stated that the specific choice of the shape function 
indeed has a significant influence for the cases of interest 
in this work. 

From a numerical perspective the velocity distribu¬ 
tion according to eqs. and ( [1Q| has to be assigned 
to individual particles entering the simulation domain. 
The angular distribution is specified by the polar angle 
0 = arccos(i?j) with [3 = I/((a^ + I) [48|, and the az¬ 


imuthal angle ip = 27TRf. Rf G [0,1] are again indepen¬ 
dent, uniformly distributed random numbers. To sample 
the energy distribution the acceptance-rejection method 
is employed in accordance with reference [22|. 


IV. RESULTS AND DISCUSSION 


A. Three Dimensional Configuration 


The geometry and the discharge parameters from sec¬ 
tion [n] specify the reference case of interest for the anal¬ 
ysis of the sputtered particle transport within the given 
MFCCP configuration. For this purpose the model de¬ 


scribed in the foregoing section has been employed. 
First, the complete three dimensional geometry is con¬ 
sidered. Figure [^presents a cross-sectional slice through 
the reactor geometry zoomed to the gap region between 
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the two electrodes as indicated in figure The simulated 
aluminum density n is plotted. The target is at the top 
at ^ = 0 cm, while the opposing electrode (substrate) is 
at z = 7.5 cm. Due to complete adsorption of the alu¬ 
minum neutrals at the surfaces, the latter essentially act 
as particle sinks. On the contrary, sputtering takes place 
at the target surface only (i.e., resputtering from the 
substrate is neglected). In longitudinal direction ^ the 
aluminum density exhibits an asymmetric profile with a 
distinct maximum on the axis of symmetry, z ^ 1.85 cm 
from the target. This is on the order of a mean free path 
which under the assumption of a Maxwellian velocity dis¬ 
tribution at T = 650 K (likely underestimating the actual 
mean free path for the non-equilibrium situation) is ap¬ 
proximately Amfp ~ ~ 1.9 cm. This 

distinct maximum is an expected feature, as previously 
reported [13 E]. While all surfaces act as particle sinks, 
only the target surface also acts as a particle source. In 
consequence, the density is minimal directly at the sur¬ 
faces (in particular the electrodes). The density at the 
target exceeds the density at the substrate. A local maxi¬ 
mum is observed in the electrode gap between. For higher 
pressures the maximum is more pronounced in the center 
and shifts towards the target as the pressure is increased. 

The most significant observation from figure is that 
the density profile in the electrode gap is laterally in¬ 
homogeneous, despite the assumption of a homogeneous 
particle flux imposed over the target surface. The radial 
dependence is indeed mainly governed by the exact ge¬ 
ometry as well as the imposed fiux. More precisely, it is 
governed by the initial distribution of sputtered particles 
(that is, the spread of the polar emission angle which 
stems from the cosine distribution) as well as the lat¬ 
eral distribution of the imposed fiux (viz., homogeneous 
here). More importantly, however, it is influenced by 
the chamber volume outside the electrode gap. Its role 
as a particle sink or source changes depending on the 
degree of backscattering and the distance to the out¬ 
side (absorbing) walls: For low pressures few particles 
are backscattered and the volume acts predominantly 
as a particle sink; for higher pressures with increasing 
backscattering the volume to some degree is a particle 
source. A substantial flux of particles may reenter the 
gap region. The problem can be fundamentally viewed 
as follows: In principle particles have a chance to be lost 
to any wall (particle sink) inside or outside the electrode 
gap region. Yet, particles that have already traversed to 
the outside region also have a chance of being inwardly 
scattered. Even though it is considerably smaller than 
the chance of traversing to any outside wall it remains 
finite. This effect is naturally amplified by the nearly 
circular geometry. The fiux of particles leaving (enter¬ 
ing) the electrode gap region through the girthed area 
^girthed = ttWD 1178 cm^ is proportional to the in¬ 
side (outside) particle density. As the net flux of particles 
is outwardly directed, a gradient in the density inherently 
establishes. Consequently, the density significantly drops 
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Figure 7. Aluminum density n (above), mean radial veloc¬ 
ity Ur (center), and mean longitudinal velocity Uz (below) 
plotted over the electrode gap z for different radial positions 
R = 0,12.5 and 24 cm (solid black, dashed red and dotted 
blue line, respectively). 


in radial direction. 

The observed density gradient from the center towards 
the electrode edge is illustrated by the density and mean 
velocity profiles plotted over the electrode inter-spacing 
for different radial positions as given in figurej^ For three 
different radial positions (as indicated in figure]^ by ver¬ 
tical dashed lines) one can clearly observe a drop of the 
aluminum density towards the electrode edge. This is ac¬ 
companied by a substantial increase in the overall mean 
flow velocity. The density is relatively homogeneous for 
R < 13 cm, which is in good agreement with the experi¬ 
mental observation of a homogeneous deposition rate on 
the substrate for radii R < 12 cm [5]. The drop in density 
towards the edge of the electrode ranges from a factor of 
1.5 to 2 (between the center and the edges) depending 
on the distance from the target. As argued earlier, this 
gradient manifests in a substantial radial particle flow. 
While the longitudinal (i.e., u^) component remains al¬ 
most constant, there is a substantial increase in ur. By 
symmetry, the net flow of sputtered particles is directed 
out of the electrode gap region. Note that the flux of 
particles is conserved (that is, no particle is lost on its 
way). 

The lateral distribution of the particle flux onto the 
substrate, respectively onto the target is another measure 
for the characterization of the afore mentioned geomet¬ 
ric influence. Note that what is shown in figure is the 
fiux of particles lost to the surfaces only. It is not super¬ 
imposed with the source flux originating from the target. 
Due to the azimuthal symmetry of the reactor chamber 
being disturbed only outside the narrow electrode gap re¬ 
gion, the fiux distribution at the substrate is also nearly 
azimuthally symmetric. (The symmetry is in addition 
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Figure 8. Lateral distribution of the absorbed fluxes of alu¬ 
minum onto a) the target and b) the substrate electrode. The 
color range is from zero to a) 2.953 x 10^^ cm“^s“^ and b) 
1.214 X 10^^ cm“2s“\ 

slightly disturbed by a number of small gas inlets em¬ 
bedded into the target.) The radial profile essentially 
projects what was observed for the density. Moreover, 
now the amount of sputtered as well as deposited mate¬ 
rial can be quantified as follows: The imposed source flux 
at the target integrated over the area is TV = fj^dAFo = 
1.6 X 10^^ cm“^s“^ X 1963 cm^ = 3.14 x 10^^ s“^. On 
the other hand, the number of particles absorbed by the 
substrate electrode is TVs = 1.95 x 10^^ s“^, while for the 
target itself it is TVt = 4.61 x 10^^ s“^. Consequently, 
a fraction of 62.1 % of the original particle fiux reaches 
the substrate, while a fraction of 14.7 % is reflected back 
and reattaches to the target, and 23.2 % is lost to the 
chamber walls. The latter can be easily verified by com¬ 
parison with the effective particle loss areas. The to¬ 
tal surface area of the electrode gap region is ^gap = 
2^target + ^girthed ~ 5104 cm^. The relative contribution 
of the girthed area then is ^girthed/^gap ~ 23 %. This 
number is in remarkable agreement with the fiux fraction 
lost to the chamber walls, although this consideration 
does not account for the directionality of the sputtered 
particles as well as the additional influx of backscattered 
particles. The former manifests in the imbalance between 
the number of particles reaching the target and substrate 
surfaces, respectively (i.e., it is not ^targetMgap ~ 38 % 
each). 

The observed deposition probabilities are moreover 
consistent with experimental observations of Rossnagel 
[49] who investigated the deposition probability in a 
cylindrically symmetric magnetron discharge under com¬ 
parable conditions. For an aluminum target it was found 
that with a target-to-substrate distance of D = 5 cm and 
an argon pressure of p = 0.7 Pa (thus pD = 3.5 cm Pa is 
comparable to pD = 3.75 cm Pa in this work) the deposi¬ 
tion probability was 60 % for the opposing sample plane, 
12 % for the magnetron plane and 10 % for the side areas. 
The first two values are in excellent agreement consider¬ 
ing the given differences. On the other hand, the discrep¬ 
ancy in the side wall contribution is simply justified by 
the fact that in the work of Rossnagel the side walls do 
not comprise of all of the remaining chamber walls (as 
compared to this work). The proposed results are more¬ 



Figure 9. Line-of-sight averaged aluminum density n obtained 
from figure [^ plotted as a function of the distance ; 2 ;. 


over in reasonable agreement with the work of Turner et 
al. [50] who determined the deposited fraction of sput¬ 
tered silicon and sodium ranging between 60 % and 75 % 
(the fraction for aluminum can only be roughly estimated 
based on the respective binding energy and atomic mass). 

The longitudinal profiles of the density and mean ve¬ 
locity along the gap (in 2 : direction) are a different matter. 
Over all radii the longitudinal mean velocity component 
(i^z) increases slightly by a factor of about 1.4 from the 
target to the substrate electrode, while the radial mean 
velocity component is nearly constant (unless very close 
to the edge of the target circumference). This indicates 
two aspects: Firstly, as previously reasoned, the geome¬ 
try has a large influence on the detailed behavior of the 
particle fluxes and therefore the deposition itself (in par¬ 
ticular laterally). Secondly, the drop in density along the 
^-axis is mainly governed by a one dimensional transport 
process where particles interact with the background gas 
(i.e., scattering and thus slowing down of particles). The 
first aspect is important for reliable predictions of the 
sputtering and deposition as concerns specific details of 
the geometry. On the other hand, the second aspect is 
equally important as it indicates the principle mechanism 
of the particle transport from the target to the substrate. 

The hypothesis of a mainly one dimensional trans¬ 
port mechanism can also be argued on the basis of the 
laterally averaged aluminum density profile. By inte¬ 
grating the cross-sectional density data displayed in fig¬ 
ure along the x-axis (cf.. Fig. the line-averaged 
aluminum density profile depicted in figure (plotted 
as a function of the gap distance z) is obtained. Of 
course due to the averaging all radial information is lost. 
It is seen that the profile shape of the averaged den¬ 
sity is governed by the homogeneous central region of 
the discharge. As evident from the spatially resolved 
density (cf.. Figs. and the latter varies from a 
round shape, which is maximum at some finite value 
1.5 < 2 : < 1.85 cm for R < 13 cm (with a min/max ratio 
of roughly 0.79 x 10^^ cm“^/1.15 x 10^^ cm“^ = 69 %), 
to a nearly linear trend at R = 24 cm. Although the 
averaged profile is slightly smoothed the principle curva¬ 
ture remains round. The governing transport mechanism 
is unaltered by the averaging. The averaged profile has 
its maximum somewhere in the electrode gap (depending 
on pressure; here at 2 : ~ 1.39 cm, roughly 25 % less than 
non-averaged). Of interest in this context is the relative 
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change in values. The averaged density ranges from a 
maximum value of n ~ 1.03 x 10^^ cm“^ at z ^ 1.39 cm 
down to n ~ 0.64 x 10^^ cm“^ at the substrate elec¬ 
trode. This is a relative drop to about 62 % the max¬ 
imum value (slightly less than non-averaged). On the 
other hand, towards the target it only drops down to 
n ~ 0.98 X 10^^ cm“^, approximately 95 % the max¬ 
imum value. Note that the obtained averaged density 
profile cannot straight forwardly be compared to exper¬ 
imental data from optical emission spectroscopy (OES) 
measurements as these stem from a line-of-sight integra¬ 
tion over the optical acceptance cone of the apparatus [3] . 
In order to compare our spatially resolved model results 
with the experimental observations, the calculated data 
would have to undergo a similar analysis. A more de¬ 
tailed comparison between experiments and simulations 
focusing on the radial dependence and the three dimen¬ 
sional features of the transport of sputtered particles will 
be the topic of future work. 


B. One Dimensional Analysis 


As previously reasoned, the governing transport mecha¬ 
nisms inside the MFC CP are of mainly one dimensional 
nature (c.f., section IVA). To underpin the previously 
said, the analysis can be significantly facilitated by mak¬ 
ing use of a simplified one dimensional transport model. 
The previously applied model can be straightforwardly 
employed. Notably, the dimensional reduction leads to 
a welcomed improvement in performance (runtime of 
hours instead of days). The major achievement, how¬ 
ever, lies within the much simplified analysis. Individual 
particles are now described in a continuous ID-3V phase 
space (i.e., three velocity components are maintained but 
their position is restricted to the z direction along the 
gap). The density and mean velocity, therefore, depend 
on ^ only. The vectorial net fiux is solely longitudinal 
(T = nUzCz). As verified later, this is a consequence of 
the azimuthal symmetry of the velocity distribution of 
the Lagrangian particles. 

Analogous to the previous analysis, the aluminum den¬ 
sity as well as the longitudinal mean velocity (now equal 
to the magnitude of the overall mean velocity) are given 
in figure for the one dimensional situation. More pre¬ 
cisely what is shown are results obtained for the ref¬ 
erence case of section with T = 650 K (solid black 
line) and for a variation of the background temperature 
T = 650 K ± 100 K (twice the confidence interval of the 
temperature measurement of Bienholz et al. 0 ). A pres¬ 
sure of p = 0.5 Pa is maintained, so the corresponding 
background densities are nbg = {6.585, 5.572, 4.829} x 
10^^ cm“^. The reduction (increase) of the temperature 
has two effects: On the one hand, the mean free path 
is decreased (increased), on the other hand thermaliza- 
tion strives towards a lower (higher) equilibrium value. 
For a reduced temperature of T = 550 K (dashed red 
lines) and a correspondingly larger background density 
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Figure 10. Aluminum density n (top) and mean longitudi¬ 
nal velocity Uz (below) obtained for the one dimensional case 
plotted as a function of z. The background temperature is 
varied: T = 550 K (dashed red), T = 650 K (solid blacM 
and T = 750 K (dotted blue). The 3D result from Fig. 
(i^ = 0 cm) is given for reference (solid gray). 


the slowing down of energetic sputtered particles is more 
efficient. Consequently, the mean velocity reduces by 
roughly —10 %. The aluminum density increases (by the 
same amount) due to flux conservation and thus accumu¬ 
lation of slow particles. More importantly, the overall be¬ 
havior remains the same. That is, the density profile still 
exhibits a pronounced maximum of n ~ 1.4 x 10^^ cm“^ 
now at 2 : ~ 2.3 cm. The contrary behavior is seen when 
the collisionality is decreased by raising the temperature 
to T = 750 K (dotted blue lines). The aluminum density 
decreases by roughly —10 % because less and faster par¬ 
ticles accumulate in the volume due to fewer collisions. 
(The mean velocity increases accordingly.) Also in this 
case the generic profile of both the density and the mean 
velocity remain similar with a maximum in density at 
^ « 1.8 cm. 


As previously reported [131 HZ] 5 f he collisionality has a 
strong influence on the strength of the local maximum 
in density (for different pressures at a constant tem¬ 
perature). In the present case the relative strength of 
the maximum in aluminum density (absolute max/min) 
ranges from approximately 1.75 at T = 550 K down to 
1.38 for T = 750 K. With the same reasoning this is 
merely a consequence of the variation in collisionality 
rather than the temperature. 


It is instructive to take a closer look at the refer¬ 
ence case of T = 650 K for a detailed understanding 
of the underlying transport (Fig. 10, solid black lines). 
The principle shape of the ID aluminum density profile 
is in fact much the same as for the 3D case. A pro¬ 
nounced maximum of n ~ 1.24 x 10^^ cm“^ is observed 
approximately z ~ 2 cm from the target (comparable 
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to n ^ 1.15 X 10^^ cm“^ ai z ^ 1.85 cm for the 3D 
results on the axis of symmetry). This is again a con¬ 
sequence of a balancing of the particle fluxes originat¬ 
ing from the electrode surfaces (cf., sink/source). The 
nearly 8 % increase is from the additional contribution 
of particles which are not “lost” to the girthed area. This 
balancing is also reflected by the observance of a higher 
density n ~ 1.15 x 10^^ cm“^ at the target, compared to 
n 0.75 X 10^^ cm“^ at the substrate. It is worthwhile 
to note that the relative drop from the maximum value 
down to 60 % at the substrate electrode is quite close to 
the drop observed for the 3D situation (69 % on the axis 
of symmetry; 62 % averaged). This supports the hypoth¬ 
esis of a similarly governed transport mechanism in both 
situations. 

The line of arguments is slightly more complicated re¬ 
garding the profile of the mean velocity. Two hypotheti¬ 
cal however short explanations can be argued: 

i) A decrease of the mean velocity from the target 
towards the substrate may be expected simply be¬ 
cause the slowing down of particles scales with the 
number of collisions on the particles’ paths (gov¬ 
erned by the mean free path). The density and 
mean velocity are balanced by the steady-state con¬ 
tinuity equation V • {nu) = 0 (i.e., the flux of par¬ 
ticles has to be conserved). An increase in density 
is predicted, which is not observed. 

ii) Based on the continuity equation and given the 
density distribution along the gap it can be argued 
that the mean velocity is fully specified by the ra¬ 
tio of the imposed flux and the aluminum density 
with u = Tq/ti. This in principle gives the correct 
profile, however, scaled by the ratio of the imposed 
to the net particle flux. 

The actual behavior is in fact the same for the ID and 
3D situation: The net flux of particles is a superposi¬ 
tion of the predefined (or imposed) flux emitted from 
the target Tq and the flux of backscattered particles 
traversing in opposite direction and reattaching to the 
target. In the ID case, this flux of backscattered par¬ 
ticles can be quantified to Tt = 3.337 x 10^^ cm“^s“^, 
which makes about 21 % of the initially imposed flux. 
Correspondingly the flux of particles that reach the sub¬ 
strate is Tg = 1.266 X 10^^ cm“^s“^, equal to a deposi¬ 
tion probability of 79 %. Expectedly both numbers are 
slightly higher compared to the 3D case, because now 
particles can either be lost to the target or the sub¬ 
strate (To = Tt + Tg). The ratio of the fluxes to the 
electrodes, however, is comparable: Tt/TgliD = 0.26 vs. 
Tt/TglsD = 0.24. Now, to explain the mean velocity rise 
two aspects are important: Firstly, the substrate strictly 
acts as a particle or flux sink. The reversely directed flux 
originates from the volume due to backscattering (i.e., it 
vanishes directly at the substrate). Secondly, the density 
as well as the mean velocity are fully specified only by 
the sum of both fluxes. As both fluxes are of opposite 
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Figure 11. Aluminum velocity distributions for the different 
velocity components and a background temperature of T = 
650 K at a pressure p = 0.5 Pa. For comparison a Maxwellian 
velocity distribution of temperature T = 1200 K is plotted 
(solid gray line). 


sign they compensate each other. When super-imposed 
they are truly subtractive, not additive. (Note that the 
density as a scalar quantity indeed behaves truly addi¬ 
tive.) This explains why the net flux as well as the mean 
velocity are overestimated from the previous argument. 

To prove this hypothesis it is instructive to investigate 
the underlying velocity distribution function (VDF) of 
the Lagrangian particles. Figure presents the three 
corresponding velocity distributions within three differ¬ 
ent regions inside the electrode gap. For comparison all 
distributions (in their respective spatial interval) are nor¬ 
malized such that = 1. Note that the in¬ 

dividual particle velocities are denoted by Uk to avoid 
confusion with the mean velocity u\^. Certain aspects are 
of importance with respect to the particle transport: 

i) The simulation model is one dimensional so by sym¬ 
metry the velocity distributions in lateral direction 
/(ux) and f{vy) are expected to be indistinguish¬ 
able (within the statistical error). Consequently, all 
but the z-component of the mean velocity Uz dis¬ 
cussed earlier actually vanish. As evident from the 
two topmost subfigures of figure pT] this is strictly 
satisfied. There is no net fiux of particles in any 
other direction than :hz. 

ii) The cosine angular distribution imposed at the tar¬ 
get boundary specifies the initial velocities with 

> 0, but non-zero and Vy. The mag¬ 
nitudes follow from the high energy tail of the 
modified Thompson distribution. Compared to a 
Maxwellian distribution the lateral VDFs f{v^) and 
f{vy) exhibit slightly overpopulated tails (only a 
fraction of particles is thermalized). Indeed the 
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VDFs of only the backscattered particles would be 
much better represented by a Maxwellian distri¬ 
bution, due to collisions with the background gas. 
This can be reasoned from the lateral VDFs for 
particles residing in different intervals within the 
electrode gap. Most of the particles that reside 
close to the target are truly original (they have un¬ 
dergone fewer collisions with the background), thus 
preserving the original shape of their VDFs. Con¬ 
sequently, the agreement with a Maxwellian distri¬ 
bution is worse: it is slightly broader (solid black 
line) compared to regions further away from the 
target (dashed red and dotted blue line). The sim¬ 
ilarity of the tail of the distributions in all three 
regions stems from the initial cosine angular distri¬ 
bution as well as the comparably small number of 
scattering events. Note that even the slow particles 
are far from equilibrium with the background gas 
with temperature T = 650 K. In the steady situa¬ 
tion their distribution can be approximated with a 
Maxwellian temperature of T = 1200 K (solid gray 
line). 

iii) The situation is different for f{vz), the distribution 
of the longitudinal velocity component. The indi¬ 
vidual particle positions are constrained to one spa¬ 
tial dimension, while three particle velocity compo¬ 
nents are maintained. Consequently, the degrees of 
freedom for the mean quantities of the system are 
restricted to this direction (i.e., the ^-direction). 
By symmetry only n and are allowed to establish 
distinct profiles along the gap distance. The latter 
simultaneously allows for a substantial variation of 
f{vz) for different regions inside the electrode gap. 
The observed distributions can be explained as fol¬ 
lows: Close to the target (solid black line) the flux 
of newly injected “original” particles is largest. At 
the same time the flux of backscattered particles is 
largest too, because these originate from the whole 
gap region behind. Consequently, in this region 
the flux of most “ideal” Thompson distributed par¬ 
ticles (positively directed) interferes with the flux 
of most-thermalized backscattered particles (nega¬ 
tively directed). Moreover, slow (scattered) par¬ 
ticles reside in a given volume for a longer pe¬ 
riod of time than fast particles and therefore seem¬ 
ingly accumulate in the volume. The described 
phenomenon can be well observed in /('T’z)- A 
substantial contribution of backscattered particles 
with Vz < 0 is observed. At the same time, the 
part of the distribution with Vz > 0 is the super¬ 
position of the “original” Thompson distribution 
(with its high energy tail) and the contribution of 
positively directed scattered particles. The latter 
contribute most to the slow particle component. 
Directly at the target surface, however, they ac¬ 
tually vanish (as there is no volume to originate 
from). A contrary effect can be observed in prox¬ 
imity to the substrate (dotted blue line). There, 


the tail of the distribution is most “distorted” (al¬ 
though it is actually quite unmodified due to the 
low collisionality). The flux of positively directed 
unaltered particles is smallest, while its scattered 
equivalent is largest. The drop in the negatively 
directed contribution comes from the fact that no 
particle originates from the substrate (a particle 
sink). Its contribution to the flux is smallest (it 
vanishes directly at the surface, assuming no re¬ 
sputtering) . This explains the increase of the mean 
velocity towards the substrate previously observed 
and not understood. It is an immediate result of 
the interference of the positively and negatively di¬ 
rected particle fluxes. The mean velocity increase 
essentially results from the particle scattering in¬ 
side the gap volume. 


V. CONCLUSIONS 

The goal of this work is to provide sound insight into the 
not well understood physics of particle transport in large 
area capacitive sputtering devices. Several novel aspects 
have been addressed: i) The no-time counter idea has 
been adopted to the test multi-particle method, ii) A 
consistent set of scattering parameters to be used with 
the Ml collision model has been obtained, iii) A modified 
angular dependence for the initial energy distribution of 
sputtered particles was proposed, iv) Finally, using the 
advocated numerical simulation the governing transport 
mechanism of sputtered particles within the MFC CP re¬ 
actor has been explained. The analysis was initially per¬ 
formed for the fully resolved three dimensional geometry. 
Therein, a strong influence of the geometry on the lateral 
transport has been observed, having immediate influence 
on the spatial distribution of coating formation on the 
substrate, the walls and the target itself. With respect to 
the longitudinal transport it has been postulated that the 
peculiar essence of the particle transport can be readily 
pictured in one spatial dimension. To justify this hypoth¬ 
esis a one dimensional simulation model was utilized. It 
has been shown that despite the drastic simplification the 
dominant transport mechanism is retained, while allow¬ 
ing for a much easier interpretation of the fundamental 
aspects. On the basis of the spatially resolved velocity 
distributions of the individual particles the observed den¬ 
sity and mean velocity profiles have been explained by 
the interference of the intrinsic fluxes of “original” and 
scattered particles. This result is particularly valuable as 
it provides an insightful interpretation of the underlying 
transport mechanism, which can be directly related to 
the three dimensional situation. An in-depth comparison 
of numerical simulations of the three dimensional situa¬ 
tion with spatially resolved experimental measurements 
will be the subject of a future work. 
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